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Abstract 

We show how to perform a resummation, to all orders in perturbation the- 
ory, of a certain class of gauge invariant diagrams in Lattice QCD. These 
diagrams are often largely responsible for lattice artifacts. Our resummation 
leads to an improved perturbative expansion. Applied to a number of cases 
of interest, this expansion yields results remarkably close to corresponding 
nonperturbative estimates. 
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I. INTRODUCTION 



Ever since the earliest days of lattice field theory, one problem present in most numerical 
simulations has been the calculation of corrections induced by renormalization on Monte 
Carlo results. Although this notorious problem has not as yet been adequately dealt with, 
several methods have been used to address it: To begin with, perturbation theory provides in 
principle a methodical means of calculating, order by order in the coupling, renormalization 
functions, operator mixing coefficients, etc. Its drawbacks lie in its asymptotic nature, 
and in that it is a formidable task on the lattice, which places severe limitations on the 
order to which it can be carried out; indeed, at present, exact calculations in perturbative 
lattice QCD reach only two loops (for 2-point diagrams) and three loops (for vacuum 
diagrams) [Q. Various nonperturbative, numerical approaches to renormalization functions 
have also been devised and there has been recent progress both in their range of applicability 
and in their precision Finally, much effort has also gone in studying improved actions 

(which may, among other advantages, show improved renormalization behaviour) and 
improved/boosted perturbation theory |10|. 

In this paper, we present an improvement of lattice pertur- 
bation theory, which results from a resummation to all orders of 
a certain class of diagrams, dubbed "cactus" diagrams. Briefly 
stated, these are tadpole diagrams which become disconnected 
if any one of their vertices is removed (see Figure 1). Our 
original motivation was the well known observation of "tad- 
pole dominance" in lattice perturbation theory (see, e.g. |jll|). 
This observation must clearly be taken with a grain of salt: 
One-sided inclusion of tadpoles can ruin desirable partial can- Figure 1: A cactus 
cellations between tadpole and non-tadpole diagrams; worse, 

their contribution is gauge dependent. The class of terms we propose to resum circumvents 
the latter objection since, as we shall see, it is gauge invariant; it also overcomes the former 
objection in known cases. 

Cactus resummation may be applied either to bare quantities or to quantities which have 
been calculated to a given order in perturbation theory; thus contributions which are not 
included in the resummation can be reintroduced in a systematic manner. 

In Section |I| we present our calculation, leading to expressions for a dressed propagator 
and dressed vertices of interest; some derivations and technical details are relegated to the 
Appendices. In Section |T|, we proceed to use these expressions to calculate various renor- 
malization functions, and compare our results with other methods: We find a remarkable 
improvement in many cases. 

II. THE CALCULATION 
A. The dressed propagator 

Consider the standard Wilson action for SU{N) lattice gauge fields: 
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^ = l^Retr(l-f/°J (1) 

Ux,^iu is the usual product of link variables around a plaquette in the ii—v plane with origin 
at X ; in standard notation it reads: 

f/° = Q^9oAx,^ ^igoA^+^^v ^-igoAx+v,^ ^-igoAx^v _ J^a 

By the Baker-Campbell-Hausdorff (BCH) formula we have: 

Ux,,^u = exp {igoi^x,^ + - A^+^^f, - A^.^) + O(^o)} 

^ exp {tgoF^^l + iglF^l, + iglF^l, + O(^7o') } (3) 

The diagrams which we propose to resum to all orders will be cactus diagrams made of 
vertices containing F^']^^, ■ Let us see how such diagrams will dress the gluon propagator; we 
write: 



where the one-particle irreducible piece is given by the recursive equation: 



o 
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+ + 
+ 



■ (5) 

Now, the fact that the vertices involved in the above contain only -FJ^^^ implies that the 
longitudinal parts of all propagators will always cancel. As we will see, this fact will lead to 
the result that the effect of dressing is the same in all covariant gauges. We will thus denote 
by a thick (thin) solid line the transverse dressed (bare) propagator. 
From Eq. (|^) there follows: 

= w{go) (6) 



Indeed, the dressed propagator will become a multiple of the bare transverse one, where the 
factor w{go) will depend on gQ and A^, but not on the momentum. Let us now turn the 
diagrammatic relations (§, |^) into an algebraic equation for w{gQ); from Eq. (|^) we have: 
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(l + ^(,o) + w{gof + ...) = —. ^-1^ (7) 



and from Eq. (^) we find: 



It is crucial to verify at this stage that all diagrams contained above appear with the same 
combinatorial factors as in the ordinary perturbative expansion; this is indeed the case. 

To proceed, we must evaluate the generic tadpole appearing in Eq. (||); this comes from 
an n-point vertex of the action, in which n — 2 lines have been pairwise contracted. Before 
contraction, the vertex reads: 



.g«(giH \-qn)x ^^^rpai'-pa2 rpa„^ ^g^ 



(g^ = 2sin(g^/2) ). At contraction there will be (n — 2)/2 loop integrations giving: 

A^ = i (10) 



For the contraction of the SU{N) generators we first define and evaluate F{n; N), which 
is the sum over all complete pairwise contractions of trlT^^T^^ . . .T""}: 

F{n; N) = ^^^7^^ E 6.,., . . . tr {t^(-) T^(-) . . . T^(-)} (11) 



(F(2n + 1; A^) = 0; Sn is the permutation group of n objects). In Appendix |^ we calculate 
the generating function of F(n; A^): 

00 n Jn 

G{z-N)^Y.^nn;N), whence: F{n- N) = —G{z; N)\^^^ (12) 

We find: 

G{z; N) = e^'(^-i)/(^^) L]v_i(-^V2) (13) 

[Lp are Laguerre polynomials). In the present case, two legs are left external, so that the 
color contraction gives: 

nF{n;N) 

Substituting (^) and (|T3|) in Eq. (^, we obtain for the tadpole: 
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J2 / <lqidq,_('ii^A'^(qi)-qi,Al{qi)){q.2^Al{q2)-q2^A 



We can now sum up all terms in Eq. (H); we obtain: 

1 2{tgor 1 



(15) 



w{go) 



V - "^"^"^ i Ffn- A^V 



(16) 



n=4,6,8. 



F(n+ 1; AT) 



2(«fl'o 



[2-2M(7o)]'/' + l 



Comparing with the definition of G{z;N), Eq. (0), we see that the expression in curly 
brackets above is simply G'{z; N), the derivative of G{z; N). Eq. (|16D now reads: 



From our result for G{z; N), Eq. (p!3|), we see that: 
zG'(z;iV) = e^'(^-i)/(^^) 



(17) 



N-1 ^ Z^ 2-2 



(18) 



This allows us to make explicit Eq. (|T^): 



u e 



-u{N-l)/{2N) 



N 



L],_^{u) + 2L 



2 I 

Af-2' 



U 



9l (iV^-1) 



(19) 



u{go) 



9l 



4(l-u<(7o)) 



Given (^o, A^, this equation can be solved numerically for u[gQ) and, subsequently, for 
w{go). The region in go for which a solution exists contains the whole range of physical 
interest. Indeed one finds a solution in the region 

< c/o < 16/36^/2 ^ 3 23 for = 2, and 
< (^o — 1-558 for A^ = 3. In Figures 1 and 2 we plot the left hand side of Eq. (|1^) versus 
u, for SU{2) and SU{3) respectively. 



B. Vertices from the action 

The 3-point vertex of the action can be dressed to all orders in a manner similar to 
Eq. (D. We have: 




+ 



1 



1 - w{9o) 



1 



[l-wigoW 



+ ■ 



(20) 



The calculation is described in Appendix O. The result turns out to be very simple: 
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where ^(5^0) is the quantity calculated previously. 

Vertices with more lines can be treated similarly; however, the dressed vertex in these 
cases is not merely a multiple of the bare one, which tends to complicate matters somewhat. 
Since we will not need such vertices for the numerical results of Section |T^, we only present 
some relevant formulae in Appendic 0. 



C. Other operators 

Various lattice operators can be dressed by the same procedure. Let us take as an 
example a typical operator involving gluons: 

O = Etr{f/if/2 . . . f/4 = E {e'^"^} (22) 

where the sum runs over the Lorentz indices involved. Using the first order BCH expansion 
for (5, we can write once again for the 2-point tadpole built out of an n-point vertex: 

A{tg,Y-' nF{n-N) 




(23) 



Here, a is the value of the one-loop momentum integration coming from the contraction of 
Q with itself; it is a pure number which depends on the operator under consideration. For 
example: 

O = S, a = as before (24) 

O = E e^'""'HU.,^.uU.,p„}, a = 1 (25) 

0=Y: HU^.AU.,up. U^,p,]}, « = ^ - Jd^^c^ 0.85332 (26) 



and so on. The complete resummation of cactus diagrams then leads to: 



e-^(A^-i)/{2iV) 



1 



x=(gf,a)/(2~2wigQ)) N^-l 



It turns out that 3-point bare and dressed vertices are related by the same proportionality 
factor as the 2-point vertices, Eq. (PTf). 

For the topological charge density operator of Eq. (p5|) an alternative resummation is 
possible by using the BCH expansion as follows: 
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Keeping the first order terms in F^^^^, and -F^.po- the complete resummation leads to the 
simple result: 

— ®— = — ®— ■ [1 - w{go)]' (29) 

The square in the above expression can be traced to the fact that the operator is composed 
of two mutually othogonal plaquettes. 

Cactus resummation can be also used to estimate the perturbative vacuum expectation 
value of an operator: 




D. Other representations 



The calculation performed above can be generalized to encompass several other cases, 
e.g. operators involving higher representations for gluons. To illustrate this, we consider the 



class of variant actions proposed some time ago 



2 ^ V "''"'7 ' 2 ^ V iV2-l 



tr^f/x, 



(32) 



Here /3 and Pa are adjustable parameters and t^AUx^^iu denotes the trace of a product of 
links in the adjoint representation, around a plaquette. 

The calculation proceeds as before; the one new ingredient we need is FxA]{n] N), the 
sum over all complete pairwise contractions of trjT'^^T'^^ . . .T"^"} (T°: generators in the 
adjoint representation). In Appendix we compute GxA){,z; N), the generating function for 
FAdj • In terms of G^A] , the equation for the factor Wvar(5'o) which dresses the propagator 
now becomes: 



4^^'^^'^) + ^(i^^^Ad,(^;iV) 



iV2 



(33) 



[2 - 2w,U9o)V^' 



where: 



9^0 



/3 PaN 
2N iV2-l 



(34) 
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It is straightforward to solve Eq. (^) numerically for Wvar(fl'o)- 

We conclude this Section by noting that the extension to vertices with fermions is im- 
mediate. First of all, vertices coming from the Wilson fermionic action stay unchanged, 
since their definition contains no plaquettes on which to apply the linear BCH formula. We 
will see how this affects corresponding renormalization functions in the next Section. For 
more complicated fermionic vertices, such as those of the clover action, cactus resummation 
proceeds in precisely the same manner as Eqs. (ESI, B^. 



III. SOME APPLICATIONS 

In this section we apply the resummation of the cactus diagrams derived for the Wilson 
action to the calculation of the renormalization of some lattice operators. Approximate ex- 
pressions for these renormalizations on the lattice are obtained by dressing the corresponding 
one-loop results. We will consider here operators whose anomalous dimensions are zero. A 
consistent, as well as physically motivated, means of implementing the cactus dressing is to 
apply it to the one-loop difference between lattice and continuum contributions that deter- 
mine the renormalization, and not only to the lattice part. Cases with nonzero anomalous 
dimension can be dealt with in an analogous manner, by setting the scale /i = 1/a and 
dressing the finite renormalization coefficients as before. 

As a first example we consider the calculation of the lattice renormalization Zlg^) of the 
topological charge density operator 

Q(^)=- n,jon 2 E ^"^^Mf/.,..f/.,pJ- (35) 

Z{g'^) is a finite function of Qq] it approaches one in the limit go ^ 0, and is much smaller 
than one in the region go — I, where Monte Carlo simulations using the Wilson action are 



actually performed. A non-perturbative numerical calculation using the heating method [13 
has produced the estimate |1T^ , [I5[| 



Z{g^ = 1) = 0.19(1) for SU{3). (36) 

In this case few terms of the perturbation theory in gQ can hardly provide an acceptable 
estimate of Zlg^) for go — 1 without some kind of resummation. 

In perturbation theory Z^g^) has been computed to O^g^) |16| with the result 

Zigl) = l + z,gl + Oig',), (37) 
zi = N - - - At - 0.15493 

Numerically zi ~ —0.908 for SU{3), and zi ~ —0.536 for SU{2). So perturbation theory 
to 0((7q) would give Zi^gl) ~ 0.092 for SU{?>), which is very far from its actual value, 
Eq. (^). In order to obtain a better approximation, we perform a cactus dressing of the 
one-loop calculation. The tree order is dressed by mere use of Equation (p9|) . One can now 
dress the one-loop contributions (for the details of the standard perturbative calculation see 
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Refs. |16| , p!7| ). Using Eqs. (0), (pi]) and p9|), and a simple combinatorial counting applied 
to the diagrams contributing to Z{g'^), one arrives at the following expression 

Z^gl) ^ [1 - w{go)f + [1 - w{go)] (^zi + 9l (38) 

The quantity (2A^^ — 3)/12N must be added to zi to avoid double counting, since such 
a contribution is already incorporated in the dressed tree-order approximation. Solving 
Eq. (|rup for = 3 and go = 1 one finds 

1 - ^'(^o = 1) = 0.749775. (39) 

Thus from Eq. ( p8|) one obtains Z{gQ = 1) ~ 0.193, which compares very well with the 
numerical result (^6]) . Further confirmation of the validity of the approximation ( |38D comes 
from a comparison with available data for SU{2) in the range 2.45 < (5 < 3.0 {(3 = ^/ g^) 
obtained by the heating method [|l8l, as shown in Table |. The agreement is remarkable. 



We wish to point out that other improvement recipes, such as those proposed in Ref. [10 
consisting in a redefinition of the bare-coupling, do not help in this case. For example one 
recipe entails use of 

r - T^y (40) 

(where [/° is the plaquette) as the expansion parameter for = 3. In many cases this recipe 
represents an improvement. However, substituting the value of g^ corresponding to (^q = 1, 
i.e. g^ ~ 1.68, in Eq. (p^), one would obtain Z{gQ = 1) ~ —0.54, that is much worse than 
the plain one-loop approximation. Similarly, a change of coupling constant and momentum 
scale, in the manner of Lepage and Mackenzie |T0|, also leads to a wider discrepancy in this 
case: indeed, the corresponding value of a{q*) (defined in |jlO|), turns out to be too large. 

One can also apply cactus resummation to the lattice renormalization of fermionic op- 
erators. Let us consider the local non-singlet vector and axial currents = V'A"7^?/' and 

= ijjX°''jfj,'y5ip. The lattice renormalizations of these operators, Zvigo) ^Aido) respec- 
tively, are again finite functions of go- In perturbation theory and for SU{3) one has |]Ty| 

Zv,A = l + zv,A9o + 0{gt), (41) 

where zy — —0.17 and za — —0.13. Thus, at g^ = 1 one-loop perturbation theory gives 
Zv{gl = 1) ~ 0.83 and ^a(5'o = 1) — 0.87. For these fermionic operators, one may use 
cactus resummation to dress the gluon propagators appearing in the diagrams contributing 
to one-loop order, according to Eq. (|^). This procedure leads to 

Zv,A-l + Zv.A , y (42) 

1 - w{go) 

At g^ = 1, this gives: Zy ~ 0.77 and Zj^ ~ 0.83. One may compare these numbers with 
those obtained in nonperturbative calculations based on the Ward identities 0. The only 
limitation of this method is due to scaling corrections, which turn out to be rather large at 
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(7o — 1 in the case of the Wilson lattice formulation. Depending on the matrix element one 
looks at, at (yfg = 1 one finds values ranging from 0.57 to 0.79 for Zy and from 0.72 to 0.85 for 
Zp^ [p0|-|2^ (see also Ref. I^S] for a review of these results). Other methods of improvement 



see, e.g., also p3[ for a partial review), using various boosting procedures, result in 



numbers ranging from 0.63 to 0.71 for Zy and from 0.72 to 0.77 for Za- Hence, a conclusive 
comparison is not possible in these cases. A better numerical situation occurs when one 
considers the clover action for which scaling corrections are largely reduced in the region 
where Monte Carlo simulations are performed and precise estimates can be obtained using 
the Ward identities (see Ref. |]23| and references therein). An application of our cactus 
resummation to the clover action would require the dressing of the new fermion-gluon three- 
point vertex. This point is under investigation. 

In conclusion, the above examples show that the resummation of cactus diagrams leads 
to a general improvement in the evaluation of lattice renormalizations based on perturbation 
theory. A combination of this method with improved actions is expected to give a reliable 
evaluation of renormalization functions, which can complement corresponding nonperturba- 
tive estimates. We hope to return to this issue in a future publication. 



ACKNOWLEDGMENTS 

It is a pleasure to thank Paolo Rossi for very useful discussions. H. P. would like to 
acknowledge the warm hospitality extended to him by the Theory Group in Pisa during 
various stages of this work. 



APPENDIX A: CALCULATION OF G(z; iV) 

We wish to calculate F{n;N), the sum over all complete pairwise contractions of 
tr{T"iT'^2 _ _ _ 2-a„|_ Yor even n, F{n; N) is defined by: 

N) = ^^^7^^ _E ^a,a, . . . tr {t^(-) T^(-) . . . T^(-)} (Al) 



PeSn 



{F{2n + 1; A^) =0). Sn is the permutation group of n objects, and are an orthonormal 
basis for su{N) in the fundamental representation, tr{T"T^} = 
We define: 

M = e^T'', = (a = 1,...,A^^-1), reIR (A2) 

Then F(n; A^) can be written as: 

F{n; N) = ^ / R tr{M"} = ^J [dM] e"^'' tr{M"}, ^ = J [dM] e"*"' 

(A3) 

(The normalization Af is redefined below as necessary, to ensure that -F(0; A^) = A^ remains 
valid). By definition, [dM] = HadO"' is the integration measure over traceless Hermitian 
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matrices. When the integrand is invariant under similarity transformations, as is our case, 
"angular" integrations can be performed, leaving behind an integral over the N eigenvalues 
A, 121: 



■'^ i i<j i \ i / 



At this stage, it is convenient to introduce the generating function for F{n; N): 



oo 



AT) ^ 5] -F(n; AT) 



n=0 



n! 



By Eq. 



we have: 



i ■' •'^ k<l rn 



N 



k<l 



To simplify the exponents we shift the A's so as to keep their sum equal to zero: 



A' - A + ^ 



a: = a,+;;^ 



z 

2N 



(A5) 



(A6) 



(A7) 



Giz;N) = Nf^[ n (A.-A.f][n(Ai + f-A,.r]5(EA.)e-i:.^?-^(--)/(^-) 

(A8) 

The (5-function can now be easily eliminated, using the exponential representation 
/ da exp(za J2i Aj): 



G(z;iV) = iVe^^(^-i)/(^^)/^^^^[ n (A.-A,)^] [ H (Ai+f -A 



^zHN-1)/{4N) 



n(A.-A,f][n(A.+|-A 

-'^ i^fc<« k^i ^ 



g-^^(A,-*a/2)2-a2 7V/2 

e-S.^? (A9) 



Let us isolate the integral over Aj (i 7^ 1): 



N) = N e^^(^-i)/(4^) | dAi e"^? e(^^+^/2)' [ | ( g (A',- A;)2)e^ 



(AlO) 



The integral in square brackets, involving the Vandermonde determinant ^k<i{^k ~ A;), 
equals 
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N-1 



d \ 



We thus obtain: 



,(Ai+^/2)2 



J! 



1 2 



h '\ 2j 

in terms of the Hermite {Hj) and Laguerre (L^) polynomials. 



(All) 



n 2 



-(Ai+z/2)2 



(A12) 



APPENDIX B: PROOF OF EQ. (|2l]) 



To prove Eq. pl|), we must first evaluate the j-loop tadpole diagrams appearing in 
Eq. (pO[). Contracted legs come from F^^^^, while external legs necessarily originate from 
^x]j.u ^^"^ -^if/li/ • "^^^ corresponding vertex comes from 



^ ^ tr (1 - exp(25(oF^,^^)) 



(Bl) 



taking 2j + 2 powers from the exponent. We have: 




(2j + 2)! 



2i + 2 

Ar2-1 



F(2j + 2;iV) 



2, 



(B2) 



The first factor above is the ratio of Taylor coefficients for the vertices on the l.h.s. and 
r.h.s.; the factor in square brackets is the outcome of color contractions; the factor (1/2)-^ is 
the outcome of j 1-loop integrations. Combining Eqs. ( P2|) and (^) we find: 




2,^! [2-2wigo)]^/^ dz 




■G'i 



2[2-2w{go)Y/^ 



[2-2w{go)Yr 



(^^7o)(iV^-l) 




[1 - w{go)] 



(B3) 



In the last equality, use was made of Eq. (IT 
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APPENDIX C: DRESSING THE 4-POINT VERTEX 



The bare 4-point vertex contains parts coming from trlFj^^^^Fj^^^^Fj^iz-^i^i^}, 
trjFj^^^Fj^^^}, trjFJi'^j^Fj^^j^}. In general, these are expected to dress differently, thus yield- 
ing a dressed vertex which is not merely proportional to the bare one. In all other respects, 
this calculation is a direct generalization of the 3-point vertex case. 

We will not present the final expression for the dressed 4-point vertex, since we will not 
be needing it in Section |ITTt rather we evaluate the one new ingredient present in this case: 
The sum over all pairwise contractions of tr with 4 legs left external. The result 

can be written as: 



('ix(n-4)/2 p(l)ap(l)6p(l)cp(l)d 

^2^' X,^lU X,flU X,fJLV X,fJ,U 



abed 



(CI) 



where T'^*^'' is necessarily of the form: 



T 



•abed 



a ■ ii{T''{T'T^T'' + permutations)} + /5 ■ {5'''5^'' + 5'"' 5''' + 5""^ 5^^^ 



(C2) 



We must compute a and [3 for generic n, A^. 

The tensors multiplying a and (3 above are in general independent, except for the cases 
= 2, = 3. One way to see this is by taking the scalar product of T"-^'^'^ (a real tensor) 
with itself: 



nabedrpabed ^ ^(^p^2 _ ^-^^^^2 

= 3(iV2-l)(A^2 + l) 



+ ap — — — - + 



{[3 + 



8Ar2(iV2 + l) 
2A^2 _3 



N{N^ + 1) 



,(iV^-9)(iV^-4) 



(C3) 



-la. 



The above can vanish only for N = 2 or N = 3, ii P 

To compute a and /? we further contract T"''^'^ with either -.d^^d^^: or ■.^{T'^T^T^T^}:, 
to arrive at two relations for a and /3: 



4!(Ar2 



2N 



[(2A^2-3) a + 2N{N^+1) p] = n{n-2) F{n; N) 



(C4) 



lem " + 4iV(2iV'-3) P] = ^ F{n+2- N) - ^^^^^^^^^ F(n; N) 



n(n-l)(iV2+n-3) 
8iV2 



F(n-2; A^) 



The solution of these linear equations gives the required expressions for a and p. 

For A^ = 2, A^ = 3, Eqs. (|C^ ) are linearly dependent, as expected. Since the two tensors 
in T"-^'^'^ are proportional to each other in this case, we can set a = 0. Then, from Eq. (|C4|) 
we find: 



a = 0, P 



n{n-2)F{n]N) 
4!(Ar2 - 1)(A^2 + ly 



{N = 2oi N = Z) 



(C5) 
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APPENDIX D: RESULTS FOR THE ADJOINT REPRESENTATION 

We calculate: 

^Adj(n; N) = ^^^7^^ _E ^a.a, 5a3a, • • • 5a„_,a„ tr {T^^^^^ T^^^^^ . . . T^^^)} (Dl) 



PeSn 



Here, T"^ denote su{N) generators in the adjoint representation. We can relate them to the 
fundamental representation using the standard decomposition: 

{N) ® (N*) — y (iV'-l) © (1) (D2) 

In terms of the generators, this says that there exists a unitary matrix U such that: 

f/"^ (T" ® U + U ® T"*) U = T^® (0) (D3) 

Using an A^^ x (A^^ — 1) projector P (the (A^^ — 1) x (A^^ — 1) unit matrix augmented by a row 
of zeroes): 

r" = pt[/^(T"®U + U® (T")*)f/P (D4) 
Substituting this in ( pij ) and making use of Eq. we find: 

FAdi{n;N) = ^J Ylde"" e'-^'" tr {(M ® 1 + U ® M*)"} 

= ^J llde^e-^'' J2 (^) tr{M"}tr{M"-™} (D5) 
The corresponding generating function is now given by: 

^Adj(n;Ar) = — GAdj(^;iV)|,=o 

oo ^ni+n2 1 /• -I „ 

GAdj(^;iV)= J2 /n^^"e-^^ tr{M"^}tr{M"n 

i,j '' •'^ k<l m 



A somewhat tedious integration over Am, by analogy with Eqs. (|Allj , |Al^) , gives: 



N-lN-l-n 2 1 ^2 N 

From this point on, dressing the variant action propagator proceeds just as in Eqs. (|],|15|), 
leading to 

14 



2(^^o) ro o... ^ll/2 



[2-2w^U9o)r +^ (D8) 



Substituting the definitions of G{z; N) and GAdj{z; N) in tlie above produces immediately 
Eq. (H). 
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FIGURES 




Figure 1: Plot of the l.h.s. of Eq. (|T9D versus u, for SU{2). The sohd part of the curve identifies 
the interval of Qq values for which a solution exists. 



3.116 




1.5 



Figure 2: As in Figure 1, for the case of SU{3). 
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TABLES 



TABLE L For the SU{2) lattice gauge theory we hst the estimates of Z{gQ) as obtained by 
the heating method |]l8| (h.m.), by the standard one-loop perturbative expansion (p.t.), and by 
cactus dressing, Eq. p^), of the one-loop calculation (d.p.t.). 



f3 = A/gl h.m. p.t. d.p.t. 

2I5 0.20(2) 0.125 0.219 

2.5 0.22(1) 0.142 0.233 

2.6 0.25(2) 0.175 0.259 
2.8 0.32(2) 0.234 0.305 
3.0 0.33(2) 0.285 0.347 
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